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We consider the reconstruction expected for the Fermi surface of underdoped YBa2Cu3064.a; in 
the case of a colhnear spin-density wave with a characteristic vector Q — {tt[1 ± 26], tt), assuming 
an incommensurability 5 « 0.06 similar to that found in recent neutron scattering experiments. 
A Fermi surface possibly consistent with the multiple observed quantum oscillation frequencies is 
obtained. From the low band masses expected using this model as compared with experiment, a 
uniform enhancement of the quasiparticle effective mass over the Fermi surface by a factor of « 7 
is indicated. Further predictions of the Fermi surface topology are made, which may potentially be 
tested by experiment to indicate the relevance of this model to underdoped YBa2Cu306+i:. 



The recent discovery of magnetic quantum oscilla- 
tions in high Tc superconductors provides a unique op- 
portunity to access the Fermi surface topology and the 
properties of the quasiparticles undergoing pair forma- 
tion [TmilllliJIiHIHllS]. The reports of multiple 
small pockets of carriers [3 [71 [S] and a negative Hall co- 
efficient [3l[T0] in the underdoped systems YBa2Cu306+2: 
and YBa2Cu40g suggest a possible reconstruction of the 
Fermi surface by a density- wave order parameter. In this 
paper, we calculate the Fermi surface topology consider- 
ing the case of a coUinear spin-density wave with a mod- 
ulation vector Q = (7r,7r[l ± 2S]), where 5 « 0.06 is pro- 
vided by recent elastic [11] and inelastic [12] neutron scat- 
tering experiments on underdoped YBa2Cu306+2:- In 
this model, we use the values of S and ratio Vg/t ^ 1 
(of the spin modulation potential to the in-plane hop- 
ping parameters) which represent the case of a conven- 
tional spin-density wave, rather than stripes of the form 
proposed in Refs. [Ol [H] . 

We begin by considering the unreconsructed tight 
binding approximation 

£k — £a + 2iio[cos(afc^) -I- cos{bky)] + 2tii[cos{akx 
+bky) + cos{akx — bky)] + 2i2o[cos(2afcj,) 

+ cos{2bky)] + 2tc cos ck^ (1) 

originally proposed for YBa2Cu306+a: by Andersen et 
al. |15j . In this approximation, we choose iio — 
- 380 meV, tn/iio = - 0.32 and i2oAio = 0.16 [I^ITe]. 
while the number of holes is tuned by Eq. We intro- 
duce an inter layer hopping term (2tc cos ck^) to model the 
deep corrugation of the otherwise cylindrical Fermi sur- 
face section reported in recent studies jSHnj. We further 
neglect the ortho-II potential [17], following reports of 
quantum oscillations in non-ortho-II ordered samples [7], 
and the effect of bilayer splitting [TSj HH] • Angle-resolved 
photoemission spectroscopy experiments on surface K- 
deposition-treated ortho-II ordered YBa2Cu306.5 yield 
a Fermi surface whose location in fc-space is consistent 
with these assumptions for Eqn. ([ij for hole dopings near 
p = 0.1 [19]. We plot the corresponding Fermi surface for 
A:^ = ±7r/2c in Fig. [l]i. 



For guidance on the appropriate choice oi S, we turn 
to recent neutron scattering experiments on underdoped 
YBa2Cu306-i-a;- Magnetic field-dependent elastic scat- 
tering experiments [TT] find evidence for long range an- 
tiferromagnetism with S ~ 0.055 in samples of similar 
composition {x = 0.45) to those of the lowest compo- 
sition (x ~ 0.49) in which quantum oscillations are ob- 
served |9|. Meanwhile, inelastic scattering experiments 
performed on samples of nominally the same composition 
{x « 0.5) as those in which multiple Fermi surface pock- 
ets are reported [3 H], find incommensurate magnetic 
excitations at energies as low as 6 meV at 8 ~ 0.06 ^2] 
in the absence of an applied field. 

We therefore consider the possibility that spin-density 
wave ordering occurs near b ~ 0.06 in YBa2Cu306.5 in 
sufficiently strong magnetic fields. If it is coUinear, the 
Fermi surface reconstruction is determined by multiple 
translations of by ±nQ (each of which we denote 
Ek+nq)- As in Cr [201 HI], this leads to a hierarchy of 
gaps of order 2A„i ~ 2T/,'"/t™~^ opening at the cross- 
ing points of bands Ek-i-nQ and ek-i-n±mQ, where 14 is 
the amplitude of the spin potential and t ~ t\Q. Pro- 
vided 14 ^ these gaps become vanishingly small for 
large m, implying that there is little difference between 
the orbits obtained in an incommensurate model with an 
irrational 8 and those obtained using a commensurate 
model with a rational 8 of similar value. We therefore 
choose to consider a rational (5, since it enables a com- 
plete picture of the reconstructed Fermi surface to be 
obtained from the eigenvalues of a matrix |16j . Given 
the limited precision of inelastic neutron scattering ex- 
periments on YBa2Cu306.5 |12j . we make the convenient 
choice oi 8 ^ — 0.0625, and compute the eigenvalues 
of the 16 X 16 matrix 

/£k 14 14 ... K 14 \ 
14 Ek+Q 14 ... 14 
K 14 ek+2Q ... 
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Additional terms (14) are introduced in Eqn. Q to ac- 
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FIG. 1: a Unreconstructed Fermi surface of YBa2Cu30a+a: 
for p = 0.1 according to Eqn. ([T]). b Reconstructed Fermi 
surface, as described in tlie text. The four prominent orbits 
are depicted in color (online), c Notional Fermi surface cor- 
responding to the eigenvalues of a reduced 3x3 matrix con- 
taining only the terms in the top-left-hand-corner of Eqn.([2|, 
which produces the same orbits as in (b). 



commodate a possible charge modulation at 2Q |16j . 
which often occurs in spin-density wave systems |21j . 

Figure [l]3 shows the in-plane Fermi surface cross- 
section at kz — ±7r/2c corresponding to the 16 eigenval- 
ues of Eqn. (|2| for the case of a simple collinear spin- 
density wave in which Vc — 0. Here, Vg is adjusted 
to a value of 0.15 x tio to yield an electron pocket 
(a) of similar fc-space area to that reported experimen- 
tally [TU], while £o is adjusted to conserve the hole filling 
at p = 0.1. The full paramagnetic Brillouin zone shown 
is 16 times larger than the antiferromagnetic Brillouin 
zone. Four primary orbits (a, /3, 7 and 6) are expected 
for this reconstructed Fermi surface, and are indicated 
in color in Fig. [TJd. The smallness of the higher order 
gaps 2A„>i implies that these four orbits can be ap- 
proximately reproduced by considering the eigenvalues 
of a reduced 3x3 matrix containing only on the terms 
in the top-right-right-most corner of Eqn. These 
are shown in Fig. [ij:. The multitude of very small hole 
orbits in Fig. [T|d (resulting from the overlap of many 
bands with relative translations of mQ where m > 2) 
are too small to have filled Landau levels in magnetic 
fields of the magnitude relevant for quantum oscilla- 
tions [H m [5j |6j [7l [9] . Furthermore, the gaps between 
them (2 Am ~ 21/™ /t"*-! < 2K,) are small enough to be 
completely broken through (having magnetic breakdown 
tunneling probabilities of « 1 [23 )• 

The corresponding quantum oscillation frequencies 
Fi — (Ai/A-Bz)h/eab for each of the four primary or- 
bits are plotted as a function of p in Fig. |2]i (where 5 
and Vs are fixed) to allow for possible discrepances in 
hole doping estimates [T] [5^. The frequencies are cal- 
culated from their fc-space areas Ai using the Onsager 
relation [22 , where Abz is the area of the paramagnetic 
Brillouin zone and a and b are the in-plane lattice con- 
stants. A value of ~ 8 meV is required to account for 
the deep corrugation of (the otherwise cylindrical) Fermi 
surface detected in recent experiments [5J [HI [H] , yielding 
two (maximum 'belly' and minimum 'neck') frequencies 
for each Fermi surface section. 

An appealling feature of the present Fermi surface 
model is that the smallest detectable quantum oscil- 
lation frequency predicted is the electron pocket (re- 
ferred to here as the a-pocket following the nomencla- 
ture introduced in Ref. [5]), similar to experimental find- 
ings [U m |5j [HI [9] . Furthermore, this pocket has the hght- 
est effective mass (see Fig.[2]D) and is the only one (of the 
four) not requiring magnetic breakdown for its observa- 
tion, suggesting that it will dominate both the quantum 
oscillation spectrum in amplitude [3 [8] and the Hall coef- 
ficient in intermediate magnetic fields [lOj . The remain- 
ing three prominent orbits (/?, 7 and i5 in Fig.[l]D) require 
varying degrees of magnetic breakdown to be observed, 
for which we make a rough estimate of their magnetic 
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FIG. 2: a A comparison of calculated quantum oscillation 
frequencies (lines) plotted for a range of p with those ob- 
served experimentally (circles, with error bars at 2a level). 
A finite value of ic ~ 8 meV yields two parallel lines (i.e. 
'belly' and 'neck' frequencies) for each Fermi surface section. 
Red and blue curves and symbols correspond electrons and 
holes respectively. Best agreement is obtained for p — 0.092. 
Experimental values for all the frequencies are taken from 
Ref. [3, where the wide magnetic field range enables accu- 
rate estimates, with the exception of taken from Ref. 
b Comparison of the calculated band masses nih (lines and 
left-hand-axis) with those m* reported experimentally (circles 
and right-hand-axis). A renormalization of tio from 390 meV 
to ~ 57 meV is required in order to align the calculated and 
experimental values, suggesting a possible uniform mass en- 
hancement of m* /rrih ~ 7 across the Fermi surface. Effective 
mass estimates are taken from Refs. [S1|S]. 



breakdown probabilities using 



The subscript i refers to a particular orbit of frequency Fi 
and cyclotron frequency Wc,i = eB/mi,^i, requiring tun- 
neling through a gap 2Ai to be observed [221 [26]. Of 
the three orbits (/3, 7 and 6), the hole pocket denoted 7 
requires magnetic breakdown through only a third order 
gap of order 2A3 ~ 2V^ /iIq (plus less important higher 
order gaps) , giving rise to a tunneling probability of or- 
der ~ 98 % for B « fj,oH = 50 T. The frequency 
of ~ 640 T reported by Audouard et al [5J |27] provides 
a possible candidate for this orbit (see Fig. |2|i). Since 
the 5 and [3 orbits require magnetic breakdown tunnel- 
ing through the larger first order gap 2Ai ~ their 
quantum oscillation amplitudes are expected to be sig- 
nificantly weaker than those originating from the a and 
7 orbits. The estimated magnetic breakdown probabili- 
ties are ~ 2 % and ~ 8 % respectively. Possibly 
consistent with these greatly reduced tunneling probabil- 
ities, weak features of similar frequency (w 1130 T and 
« 1650 T in Fig. |2|i) to those in the model are reported 
in Refs. [5J[S], although the harmonics of are expected 
to occur nearby in frequency [5] . 

We find that the ratio of to^ to to* in Fig. |2jD found 
experimentally is similar to that of TOb,/3 to TOb,a pre- 
dicted by the model, suggesting a possible uniform en- 
hancement of the quasiparticle effective mass over the 
Fermi surface of « 7 (relative to the bandstructure es- 
timate |15jl. On considering this renormalization factor 
over the entire Fermi surface, the total summed effec- 
tive mass of all Fermi surface sections within the an- 
tiferromagnetic Brillouin zone becomes m^-pyi ~ 7 ™e 
(where TOc is the free electron mass), compared to 
TOp « 11 Too for the single large orbit of the unrecon- 
structed paramagnetic Fermi surface. Using these esti- 
mates, we obtain electronic coefficients of the heat capac- 
ity of 7AFM = 1-46 x(to^pj^/too) = 10 mJmol~^K~^ [5] 
and 7p — 1.46 x(top/toc) = 16 mJmol^^K^^ for the re- 
constructed and unreconstructed Fermi surfaces respec- 
tively. 

One notable caveat in assessing the applicability of the 
spin-density wave model to YBa2Cu306+2; presents it- 
self in Fig. [2^ on comparing the experimental frequen- 
cies with those predicted in the model. Best agreement 
is obtained for p = 0.092 [23], which falls short of that 
p ~ 0.098 ± 0.001 estimated from the lattice parameter 
c [Il[13]. While some uncertainty in t2o/iio or in/iio in 
Eqn. ([T|) or our neglect of the ortho-II potential (which 
occurs at 8Q here) can likely account for some of the 
discrepancy, our neglect of the 2Q charge modulation 
potential may be more significant. On including Vc in 
the simulations in Fig. [3] the electron pocket shrinks or 
expands depending on whether Vc is positive or negative, 
suggesting that improved consistency with the nominal 
hole doping level could be achieved by the simultaneous 
adjustment of Vs and Ve- 
in summary, by using a value 5 ~ 0.06 for the in- 
commensurability parameter found in recent neutron 
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FIG. 3: The Fermi surfaces (shown for | of the paramagnetic 
Brillouin zone) obtained on including a significant charge 
modulation potential at 2Q, with Vc = in (a) and K = ~Vs 
in (b). 



scattering experiments, we find that the existence of a 
coUinear spin density wave in suitably strong magnetic 
fields would yield a reconstructed Fermi surface consis- 
tent with the observed multiple carrier pockets in un- 
derdoped YBa2Cu306+a; [3 ID- Four major orbits are 
predicted, with magnetic breakdown being an important 
factor for all but the smallest electron pocket. 

A key feature of this model is that the anticipated 
shape of the electron pocket is a rounded rectangle 
with an aspect ratio >w 2 and its long axis parallel to 
ky (orthogonal to the incommensurate spin modulation 
direction)— something that can be tested by dual- axis 
angle-dependent magnetic quantum oscillation or magne- 
toresistance experiments. A small electron pocket that is 
found to be circular in cross-section (or have an approxi- 
mate fourfold symmetry) in experiments would be unrec- 
oncilable with the present coUinear density wave model 
and would instead have different implications. Such 
a pocket would then suggest an evolution of 8 rather 
abruptly to a value significantly greater than 0.06 in 
strong magnetic fields [H], or a helical or spiral spin- 
density wave |29j , or an entirely different form of Fermi 
surface topology in YBa2Cu306+£!; (possibly unrelated to 
the magnetic diffraction peaks). If experiments detect 
a Fermi surface topology similar to that yielded by the 
coUinear spin density wave model considered here, indica- 
tions would be that such a form of order could be chiefiy 
responsible for Fermi surface reconstructuion in strong 
magnetic fields in underdoped YBa2Cu306+a;. Given the 
small value of that we consider in this model, a fur- 
ther implication would be that the observed quasiparticle 
mass enhancement has an origin mostly extrinsic to this 
form of coUinear density wave. 
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